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ABSTRACT 

We investigate a cross-correlation between the weak gravitational lensing field of the large- 
scale structure, k, and the secondary temperature fluctuation fleld. A, of cosmic microwave 
background (CMB) induced by Thomson scattering of CMB photons off the ionized medium 
in the mildly nonlinear structure. The cross-correlation is expected to observationally unveil 
the biasing relation between the dark matter and ionized medium distributions in the large- 
scale structure. We develop a formalism for calculating the cross-correlation function and its 
angular power spectrum based on the small angle approximation. As a result, we find that 
leading contribution to the cross-correlation comes from the secondary CMB fiuctuation field 
induced by the quadratic Doppler effects of the bulk velocity field v of ionized medium, since 
the cross-correlations with the 0{v) Doppler effect and the Ostriker-Vishniac effect of 0{v5) are 
suppressed on relevant angular scales because the linear dependence on the bulk velocity field 
leads to cancellations between positive and negative contributions among the scatterers for the 
ensemble average. The magnitude of the cross-correlation can be estimated as (kA) ~ 10~^° 
on small angular scales {I ^ 1000) under the currently favored cold dark matter model of the 
structure formation and the simplest scenario of homogeneous reionization after a given redshift 
of Zion ~ 5. Although the magnitude turns out to be small, we find several interesting aspects 
of the effect. One of them is that the density-modulated quadratic Doppler effect of 0{5eV^) 
produces the cross-correlation with the weak lensing field that has linear dependence on the 
electron density fluctuation field 6e or equivalently on the biasing relation between the dark 
matter and electron distributions. In other words, the angular power spectrum could be both 
positive or negative, depending on the positive biasing or antibiasing, respectively. Detection of 
the cross-correlation thus offers a new opportunity to observationally understand the reionization 
history of intergalactic medium connected to the dark matter clustering. 



Subject headings: cosmology:theory - cosmic microwave background - intergalactic medium - 
large-scale structure of universe 



INTRODUCTION 



The absence of the hydrogen Gunn-Peterson trough in quasar absorption spectra ( ]Gunn, fc Peterson 
1965) strongly indicates that the smoothly distributed hydrogen in the intergalactic mcdimn (IGM) is already 
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highly ionized by z = 6 (e.g., Lanzetta, Wolfe, & Turnashek 1995; Fan, et al. 2000). Recently, Becker et al. 
(2001) reported the first detection of the Gunn-Peterson test through in a z = 6.28 QSO spectrum. However, 
even a small fraction of neutral hydrogen in the IGM would result in an undetectable flux because of the 
large cross section to Lya photons, and thus further careful investigations will be necessary. These results 
also imply the presence of very luminous ionizing sources such as QSOs and high mass stars at high redshifts 
to cause photoionizations of hydrogen by z ^ 6. On the other hand, it is known that earlier reionization 
would have brought the last scattering surface to lower redshift, smoothing the intrinsic anisotropics in 
cosmic microwave background (CMB) radiation generated at the recombination epoch z ^ 1000 ( Bond et 
al. 1991 ; Sugiyama, Silk, fc Vittorio 1993[) . The recent measurements of the CMB anisotropics on subdegree 
angular scales ( ie Bernardis, et al. 2000 ; Hanany et al. 2000 ; Netterfield, et al. 2001 ) can be used to set 
an upper limit on reionization redshift such as Zion ~ 30 ( Tcgmark, fc Zaldarriaga 2000| ). Thus, under the 
currently favored adiabatic cold dark matter (CDM) scenario for the structure formation, reionization of 
the universe is expected to occur in the range of 6 ^ Zjon ~ 30. Revealing how reionization proceeds in the 
inhomogeneous universe stands out as one of the most important problems in the observational cosmology 



(e.g., see Barkana, fc Loeb 2001 for a recent review 



Secondary CMB anisotropics should be induced by scattering of CMB photons off electrons in the ionized 
medium undergoing bulk motion in the large-scale structure. It is therefore expected that the secondary 
anisotropics can provide a laboratory for studying structure formation in the reionized epoch. Kaiser (1984) 
investigated the CMB distortion induced by the bulk velocity to pure linear order of 0(v) and found that the 
induced anisotropies are small at arcminite scales resulting from cancellations between positive and negative 
shifts along the line of sight, since the gravitationally induced bulk velocities are irrotational. Ostriker fc 
Vishniac (1986) then pointed out that the coupling of the Doppler effect with spatial variations in the free 
electron density produces [iK anisotropies on arcminite scales at the order of 0(bv)^ which is the so-called 
Ostriker- Vishniac (OV) effect, because the effect is a second-order contribution of perturbations that does 



not suffer from the cancellation (see also Vishniac 1987; Efstathiou 1988; Hu, fc White 1996; Jaffe, fc 



Kamionkowski 1998). For the same reason, spatial inhomogeneities of ionization fraction associated with a 



realistic patchy reionization scenario should lead to the secondary CMB anisotropies by the Doppler effect 
(Aghanim et al. 1996; pruzinov fc Hu 199^ ; Knox, Scoccimarro, fc Dodclson 199^; Benson et al. 2001; 



Gnedin, fc Jaffe 2001). In particular, Benson et al. (2001) investigated the secondary CMB fluctuations 
under a realistic model of the reionization of IGM caused by starts in high-redshift galaxies, based on 
a semi-analytic model of galaxy formation and a high-resolution N-body simulation of the dark matter 
clustering. They showed that measurements of CMB anisotropies could actually put stringent constraints 
on the reionization processes. On the other hand, from theoretical aspects, it is not so evident to answer the 
problem which second-order effect of perturbations produces dominant contribution to the secondary CMB 
fluctuations. Hu, Scott, fc Silk (1994 hereafter HSS; see also Dodelson fc Jubas 1995) derived a complete 
Boltzmann equation including the quadratic Doppler effects O(w^) of bulk velocity that governs redshift 
evolution of the CMB temperature fluctuation field. Consequently, they found that, although the quadratic 
Doppler effect should induce spectral distortion of CMB like the thermal Sunyaev-Zel'dovich (SZ) effect 



(Zel'dovich, fc Sunyaev 1969) as a new qualitative feature, the magnitude is always smaller than that of the 



OV effect for the CDM models. 

A great uncertainty in models for structure formation in the reionized epoch is the extent and nature of 
reionization, because this requires the understanding of various complicated astrophysical processes such as 
the nature of first ionizing sources, the fraction of ionizing photons that escape the sources to become available 
for photoionization of the hydrogen in IGM, and the radiative balance between the escaped emissivity and 
the recombination rate in the inhomogeneous IGM. For example, even a naive question which high-density 
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or low-density region in the large-scale structure becomes first photoionizcd bas not still been answered 
definitely. Assuming ab initio that the distribution of the ionization fraction is given as a function of gas 
density, Miralde-Escude, Hanehnelt, & Rees (2000; hereafter MHR) suggested that the UV radiation may 
escape into and ionize the low-density region, where the recombination rates are the lowest, even though it 
is expected that the ionizing sources could be associated with the peaks in the density field. This conclusion 
has also been confirmed by the recent numerical simulations taking into account the radiative transfer 



approximately (Gncdin 2000) 



While the relation between the dark matter and ionized medium distributions is thus a crucial clue 
to understanding the reionization history connected with the gravitational instability scenario of structure 
formation, the relation seems difficult to be directly observed by methods so far proposed. Therefore, 
the purpose of this paper is to investigate the cross-correlation between the dark matter distribution and 
the secondary CMB temperature fluctuation field induced by the ionized medium in the mildly nonlinear 
structure. It is now recognized that the weak gravitational lensing effect due to the large-scale structure 



can be a unique tool to observationally reconstruct the dark matter distribution (e.g., see Bartelmann, & 
Schneider 2001 for a review). There are two feasible methods often considered in the literature; one is 



using the distortion effects on distant galaxy images that have been actually detected (e.g.. Van Waerbeke 



et al. 200C), and the other is using the weak lensing effects on the CMB map (seljak 
Zaldarriaga fc Seljak 1999; 



Takada, & Futamase 200C; Takada 2001 



Takada, Komatsu. & Futamase 1999; Takada, & Futamase 2001 



Hu 2001a 



Zaldarriaga 199£ 



Schmalzing 



b). In particular, it is expected that the method of 



using the distortion effects on the CMB map developed by Seljak & Zaldarriaga (1999; also Zaldarriaga fc 



Seljak 1999 ) allows us to reconstruct the projected dark matter distribution with high accuracy. Based on 
these considerations, in this paper we assume that the weak lensing field is a priori reconstructed by those 
methods, for example, using the future CMB data from the satellite mission Planck Surveyor. Then, taking 
a cross-correlation between the weak lensing field and the CMB map including the secondary fluctuations 
would be a challenging test. Furthermore, the cross-correlation method could improve the signal-to-noise 
(S/N) ratio of the Doppler effect compared with the original S/N in the CMB map itself if the secondary 
CMB field correlates with the weak lensing field, which is analogous to, for example, the cross-correlation 



between the CMB temperature fluctuation and polarization flelds (Hu, & White 1996; see also Peiris, & 



Spergel 2000 and Van Waerbeke, Bernardeau, & Benabed 2000 for other examples). In this paper, based on 



the small angle approximation, we first develop a formalism for calculating the cross-correlation function and 
its angular power spectrum between the weak lensing field and the secondary CMB fiuctuation field including 
the Doppler effects up to the second-order of the bulk velocity field. We then estimate the magnitude of 
the cross-correlation under the currently favored CDM models and the simplest assumption of homogeneous 
reionization history after some epoch ^ion ~ 5. 

The layout of this paper is as follows. In §H we review the relevant properties and parameters of the 
adiabatic CDM scenario for structure formation. In §H we briefly review the Boltzmann equation including 
the Doppler effect up to the second-order of bulk velocity and discuss the physical implications of the effect 
following HSS. In §^ we develop a formalism for calculating the cross-correlation function and its angular 
power spectrum between the weak lensing field and the secondary CMB temperature fluctuation field induced 
by the bulk motion of the ionized medium in the large-scale structure. Results for the currently favored 
CDM models are presented in §|^. In §^ we present the discussions and conclusions. Throughout this paper, 
we employ unit of c = 1 for light speed c. 
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2. PRELIMINARIES 

2.1. FRW Model 

The expansion rate of the universe can be described by the Friedmann equation: 

1 da 



H{t) = --17 - Ho^Jn^^oa-^ + Oao, (1) 
a at 

where a{t) denotes the scale factor, Hq — lOO/i km s^^Mpc^^ is the Hubble constant, fimo and JIao denote 
the density parameters of the present-day non-relativistic matter and the contribution of the cosmological 
constant relative to the critical density Pcro = S-ffg/STrG, respectively. Here we have set a(to) = 1 at present 
for convenience. Throughout this paper, we consider only a flat universe model as supported by the recent 
high precision CMB experiments (e.g., |de Bernardis, et al. 2000 ; Nettcrfield, et al. 2001) and then we 



have fi^o = 1 ^ ^^mo- The conformal time which we will often use in the following discussions is defined by 
dr/ — a~^dt. It is convenient to introduce the comoving angular diameter distance r(ry) with respect to i] or 
equivalently redshift z, and for a flat universe model it is given by 

K^)-^o-^ = ^^, (2) 

where 770 is the present conformal time and z is given by 1 + z — 1/a. 

It is currently believed that the reionization of the universe occurred at 6 ^ Zjon ~ 30 such that the 
reionized medium is optically thin to the Thomson scattering of CMB photons, which is characterized as 
T < 1 in terms of the optical depth, r. Assuming homogeneous ionization after some reionization epoch Zion, 
r can be analytically expressed as a function of z for a flat universe model: 



2 Ti^JC r 1 

dt'neCTT = -T;-^ Vl - fimO + f^mo(l + zY " 1 , (3) 



with 

TH = n,Q<nH^^ = 0.0692(1 - Y)n^oh. (4) 

Here (Tt is the Thomson cross section, is the mean number density of free electron, is the mean 
ionization fraction, which does not depend on time here, fibo is the present-day density parameter of baryon 
in units of the critical density, th is the optical depth to Thomson scattering to the Hubble distance today 
under the assumption of full hydrogen ionization and homogeneous distribution of baryon, and Y is the 
primordial helium fraction. We fix y = 0.24 throughout this paper. 

Again in the case of homogeneous and full hydrogen ionization, the visibility function that gives the 
probability of last scattering within dr] of 77 is expressed as 

5(77) — fe^^ — n^aTae^^ — TnHoa^^e^^ , (5) 

where overdot denotes derivative with respect to the conformal time 77. 



2.2. Perturbation Theory of the Large-Scale Structure Formation 



The evolution of the large-scale structure after the decoupling epoch ^ 1000 is characterized by the 
time evolution and scale dependence of the density contrast field d{x,t) mainly caused by the dark matter 
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distribution. The gravitational instability then leads to bulk motions of structure with the velocity field 
v(x), and we expect that the ionized baryonic medium also follows the bulk motion on large scales, leading 
to the Doppler effects on CMB photons through the Conipton scattering by electrons. Since we focus our 
investigations on the secondary CMB fluctuations up to the second-order of v, we need to consider the fields 
of S and v up to the second-order that are valid in the mildly nonlinear regime such as 5 ^ 1: 

S{x) ^ S'^^^x) + S'^^\x) + ■ ■ ■ , 

v{x) = v''^\x) + v''^\x) + ■ ■ ■ . (6) 

The quantities with superscripts '(1)' and '(2)' denote the first- and second-order perturbation quantities, 
respectively. The time dependences of these fields can be expressed in terms of cosmological parameters 



based on the (Eulerian) perturbation theory (e.g., Peebles 1980 ). In the linear regime, the time dependence 



of (5'^' is given by the growing factor as S'-^Hx,t) = D{t)S^^'){x) with 

m « H{t) f jp^, (7) 
Ja(t)[Hay 

where we have neglected the radiation energy contribution, and will normalize D{to) ~ 1 at the present time 
to for convenience in the foUowings. 

The Fourier transformation of (5^^-'(a;) is expressed as 

If the linear density fluctuations are Gaussian, the statistical properties can be completely specified by the 
power spectrum P{k) defined by 

{dkSl,)^i2nfPik)S'{k~k'), (9) 

where (. . .) denotes an ensemble average over all realizations. Note that time evolution of the power spec- 
trum is given by P{k,t) = D^{t)P{k). In this paper, assuming the adiabatic CDM scenario, we employ 



the Harrison-Zel'dovich spectrum and the BBKS transfer function (Bardeen et al. 1986) with the shape 
parameter given by Sugiyama (1995) as for the shape of P{k). The free parameter is then the normalization 
of P(fc) only, which is conventionally expressed in terms of the rms mass fiuctuations of a sphere of 8/i~"'^Mpc 
radius, i.e., ag. From the continuity equation, we can obtain the Fourier components of the bulk velocity 
filed w^'' at the first order as 

Likewise, based on the perturbation theory, we can obtain the second-order Fourier components of 
density fluctuation and bulk velocity fields in terms of 5^^^ and D (e.g., see Peebles 198(]| ): 



°k-k'' 



ik „■ f(fk' 



J±^Gik',k-k')Spil^„ 



^fii)=^-i2DD li^Gik',k-k')S'/_,'6):'„ (11) 



with 

5 ^ f ki ^2 \ fci • fc2 2 / fci • fej 



7 2 \A;2 ki J kik2 7 \ fcifc2 
3 1 /fci fc2\ fci • fca 4 /fci • k 
7^2\k^^h) kik2 ^ 7 [ kik2 



F{ki,k2) = 

G(fci, fc2) = ^ + o ^ + IT -171— + n -in— ' (12) 
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where we have neglected extremely weak dependences of the functions F and G on cosmological parameters 
rimo E^nd Qxo for simplicity, and considered the scalar type perturbations only. It is worth noting that, even if 
perturbations start from Gaussian in the linear regime, the mildly nonlinear evolution leads to non-vanishing 
third-order moments such as (S^) and {Sv'^), because those quantities are essentially the fourth-order moments 
of linear perturbation quantities like (((5*^^))^). 



3. BOLTZMANN EQUATION FOR SECOND-ORDER COMPTON SCATTERING 

Here, following HSS, we briefly review the derivation of the secondary CMB temperature fluctuation 
field induced by the Doppler effect of ionized medium up to the second order O(v^) of the bulk velocity. 
Let us first write down the Boltzmann equation in terms of the CMB photon energy p and direction 7' in 
the presence of Compton scattering between the CMB photons and free electrons that could lead to energy 
changes; 

df df dx' dfdp df 

where / is the photon occupation number and C{x,p) is the collision term by the Compton scattering. HSS 
derived the collision term that leads to the Doppler effect up to the second order 0{v^) as 



Cvix,p) = ricCTT 



(14) 



where we have ignored the term by thermal electron scattering, leading to the thermal SZ effect. The first 
term in the large brackets on the r.h.s represents the Doppler effect of 0{v) that produces the OV effect 
from the coupling with the electron density fluctuation field at the order 0{v5e). The other terms give 
the quadratic Doppler effect of O(w^). It is interesting that the effect leads to energy exchanges between 
CMB photons and electrons like the thermal SZ effect, and therefore causes the Compton y-distortion. As 
shown by HSS, we can get more physical insight into this effect by considering the limit of many scattering 
blobs or equivalently the multiple scattering limit. If we can average equation (|lj) over the direction of bulk 
velocity field v, the 0{v) Doppler effect primarily cancels in this case and then we can obtain 



C4x,p)^n.ar — ^- 



dp 



(15) 



If introducing the effective temperature T^s defined by (w^) = STes/nie (me is the electron mass), this 
collision term reduces to the same form of that in the Kompaneets equation with the thermal temperature 
Te of hot electron, which gives the thermal SZ effect. This implies that many ionized blobs distributed 
in the large-scale structure can be considered like hot electrons inside one cluster, and thus those blobs 
induce the energy exchange to the CMB photons characterized by the effective temperature Toff. Under the 
CDM scenario for the structure formation, the amplitude of Compton y-distortion induced by the quadratic 
Doppler effect is below that due to the thermal SZ effect mainly caused by intracluster hot plasma. This 
is because we can typically estimate Teff lO^-'^eV for the gravitationally induced bulk flow and Tg ~ keV 
for intracluster hot plasma, and the fact of Tg ^ Tcff leads to the conclusion, even though clusters are rare 
objects in the universe. 

In this paper, we do not employ the many scattering blobs limit for the purpose of generality as done 
by HSS. Namely, we shall no longer be concerned with frequency dependence of the temperature fluctuation 
field. Integration of equation ( [T^ ) over photon energy p leads to the Boltzmann equation up to the order 
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0{v'^) to govern the time evolution of temperature fluctuation field, A(= ST/Tcmb), along the line of sight 
with direction of 7: 

A + -f'diA = n^aTa [Aq - A + 7 • v - i;^ + 7(7 • 'i;)^] , (16) 

where we have neglected the gravitational effect and the correction due to the quadrupole temperature 
fluctuations. By solving the equation above formally, we can obtain the secondary fluctuation field induced 
by the Doppler effect of ionized medium during propagations of the CMB photons from the decoupling 
^ 1000 to present: 

rvo 

A(ryo,7)=/ dTjgiri)X,irj)[l + S,{x,tj)] [j ■ v - + lif ■ v)'] , (17) 

where 77* is the conformal time at the decoupling epoch « 1000, 3(77) is the visibility function given by 
equation , and we have used the fact that the free electron density field can be expressed as 

- nwa-Hv)Xeiv)i^ + Se). (18) 

The Xe{ri) is the ionization fraction as a function of r/, de is the density fluctuation field of electron 
along the line of sight, and nno is the mean hydrogen number density at present given by nuo = 1-12 x 
10^^cm^'^r2bo^^(l — Y)- Note that, as explicitly introduced by Gnedin & Jaffe (2001), rie can also be ex- 
pressed as Ue = Xe{l + 6x)nii — XeUiii^ + 6xix, ?;)][1 + Si){x, 77)], where Sx and 6h are the fiuctuation fields 
of the ionization fraction and the baryon density, respectively, as a function of the spatial position and time. 
It is thus clear that Se is different from Sh in a general reionization history. Note that we will deal with 
throughout this paper. 



4. FORMALISM: CROSS-CORRELATION BETWEEN DOPPLER EFFECTS ON CMB 

AND WEAK LENSING FIELD 

In this section, we develop a formalism for calculating the cross-correlation function and its angular 
power spectrum between the weak lensing field and the secondary CMB anisotropics field induced by the 
Doppler effect of the ionized medium. 



4.1. Weak Lensing Field 



The weak gravitational lensing due to the large-scale structure is now recognized as a unique tool for 



directly measuring the dark matter distribution. Seljak & Zaldarriaga (1999) (also see [Zaldarriaga fc Seljak 



1999) developed a useful method for reconstructing the projected field of dark matter fiuctuations between the 



last scattering surface and present from the lensing distortions to the product of gradients of the temperature 
filed. This method could be successfully performed by the satellite mission Planck. It has been also shown 
that a new power of temperature fiuctuations generated by the weak lensing effect on arcminute scales 
below the Silk damping scale can be directly used to reconstruct the projected dark matter field, provided 
that we have data with arcminutc-scale spatial resolution and sensitivity of ~ IO/kK — arcmin ( Seljak fc| 



Zaldarriaga 2000; Zaldarriaga 2000; Hu 2001a) and we can remove other secondary signals such as the SZ 



effect by taking advantage of the specific spectral properties or by combining the CMB measurement with 
the X-ray or the SZ effect itself. Furthermore, if measurements of the CMB polarization become possible. 
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combinations of the lensing effects on the temperature fluctuation and the E- and i3-typc polarization fields 



can improve the reconstruction (Zaldarriaga fc Seljak f99g; Hu 200fb). This is because the lensing effect 



induces the i?-type polarization on the observed sky from the non-linear coupling to the primary i?-mode, 
evin if the primary CMB contains the i?-mode only as suggested by starndard inflationary scenarios, and 
there is negligible cosmological contamination of the polarization field in the arcminute regime unlike the 



temperature fluctuations (EIu 2000a) 



The projected dark matter field is called the convergence field, k, and can be expressed in the form of 
a weighted projection of the three-dimensional density fluctuation field along the line of sight; 



driW{ri^,ri)S{x,r]), 



where W{r]^, rj) is the lensing weight function defined by 



2 r(77*j 



(19) 



(20) 



On the other hand, the lensing distortion effect on distant galaxy images can be similarly used to measure 
the relatively low redshift dark matter distribution such as z ^ 2, depending on typical redshift of source 



galaxies (Bartelmann, & Schneider 2001). In this case, the lower limit in the integration ( |19[ ) becomes the 
conformal time corresponding to the redshift of source galaxies Zs- However, we should bear in mind that the 
lensing effect on galaxies is relatively irrelevant for the purpose of this paper, because this method usually 



restricts the convergence field to low redshifts (e.g. see Van Waerbeke et al. 2000 ), while the reionization 
signal in the CMB map is imprinted at moderate and high redshifts. The weak lensing effects on CMB thus 
has a great advantage of probing the dark matter distribution up to high redshift, which cannot be attained 



by any other means (see also Takada, Komatsu, & Futamase 1999; Takada, & Futamase 2001; Takada 2001). 
In this paper, we assume that the convergence field ^(7) with respect to direction 7 is a prior reconstructed 
from measurements of the lensing effects on CMB. 



4.2. Doppler Effect of 0{v) and 0{v5e 



From equation (|T^, the secondary temperature fluctuation field induced by the 0{v) Doppler effect can 
be written as 



A(7) 



dr]g{r^)X,il + 5eh' 



(21) 



The pure linear order Doppler effect, A ~ 0{v), is suppressed on small scales because the gravitational 
instability generates irrotational flow, leading to cancellations between positive and negative Doppler shifts 



along the line of sight ( Kaiser 1984 ). For this reason, leading contribution to the secondary temperature 
fluctuations due to the Doppler effect comes from second order perturbations that do not suffer from the 
cancellation like the OV effect of 0{vSe) (|Ostriker, fc Vishniac 1986|; [Vishniac 1987|; [Efstathiou 1988t [Hu, fc 



White 1996; Jaffe, & Kamionkowski 1998) and the patchy ionization effect of 0{v5x) ( Aghanim et al. 1996| ; 
Knox, Scoccimarro, & Dodelson 1995; Gruzinov & Hu 1998; Benson et al. 2001), where 6x is the spatial 
inhomogeneities of ionization fraction as discussed below equation (|l8|). 

From equations (|l^) and (^ij), a cross-correlation between the convergence held and the secondary 
temperature fluctuation field due to the Doppler effects of 0{v) and 0{vSe) depend on the moments {5v) 
and (Svde) of perturbations, respectively, where the latter moment seems to produce some contributions in 
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the quasi nonlinear regime as discussed in § [2.2| . However, it is clear that, since those moments has linear 
dependence on the bulk velocity field i), the cross-correlation is suppressed on relevant small angular scales 
because of the cancellations between positive and negative contributions from the bulk velocity field for the 
ensemble average. Note that the two-point correlation function of the OV effect or equivalently its angular 
power spectrum arises from the quadratic contribution of the velocity field as 0{v'^S^), whereas it does not 
suffer from the cancellation (Ostriker, & Vishniac 1986; Vishniac 1987; Efstathiou 1988| ). 



4.3. Quadratic Doppler Effect of 0{v'^) 

We next consider secondary temperature fluctuations induced by the quadratic Doppler effect of the 
bulk velocity. In the case of assuming the homogeneous ionization history, the secondary fluctuations are 
divided into two contributions; the Doppler effect of 0{v^) and the density-modulated effect of 0{v^Se)- 
As will be shown later, the effect of 0{v'^Se) generally produces dominant contribution on small angular 
scales I ^ 1000, where the nonlinear structures with (5^1 arc important. From equation (|l7|), the secondary 
fluctuation field of 0{v^) is given by 

A" (7) = / drjgi7j)X, [-v' + 7(7 • v)^] . (22) 

The bulk fiow in adiabatic CDM models arises mainly from the linear regime. In the ACDM model, half 
the contribution to (v^) comes from scales k <, 0.07ft-Mpc~"'^, and the density fluctuations go nonlinear at 
k > 0.15/iMpc"\ This is the reason that the secondary CMB anisotropies by the Doppler effects of 0{v^) 
and 0{v^5e) are the order of magnitude smaller than the anisotropies by the OV effect of 0{v6e) on relevant 
angular scales. However, from the fact that there is no cross-correlation between the convergence fleld and the 
OV effect, the quadratic Doppler effects can provide primary contribution to the cross-correlation between 
the convergence field and the secondary CMB temperature fiuctuations induced by the ionized medium. 

Since we are interested in the cross-correlation on small angular scales such as I ^ 100, it is useful 
to employ the small angle approximation originally developed by Bond & Efstathiou (1987). In this case, 
the unit vector 7 in equation (|2^), which denotes the direction of line of sight, can be considered to have 
component as 7 « {6x,9y, Vl — 0"^) for |02:|,|0y| ^ 1 in the small region around the North Pole on the 
celestial sphere. It is then convenient to introduce the two-dimensional vector 0, which has 6 = {9x,9y,Q). 
Within the context of the small angle approximation, the cross-correlation function between the convergence 
field K and the temperature fluctuation fleld A" separated by the angle of 9 is deflned by 

C"^^\9) ^ (.j(7jA''t7fc))7„-7.=cose « {<S,)^^\db))^e^_0^i^e^ (23) 

where 6a and 61, are the two-dimensional vectors corresponding to unit vectors 7^ and 7^ for k and A" , 
respectively. Furthermore, if using the two-dimensional Fourier transformation as, for example, k{6) = 
/(i^Z/(27r)^K^e'''^, the cross-correlation function can be expressed in terms of its angular power spectrum 
by the Hankel transform of zeroth order as 

2 fldl 2 

C"'" (9) = J—Cr M19), (24) 

where Jo{x) is the zeroth order Bessel function. The angular power spectrum C;"" can be also expressed in 
terms of the two-dimensional Fourier components of k and A" as 

{niAi:)^Cr\27TrS'{l + l'). (25) 
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From equations and (|22|), the cross-correlation function (|2^) can be further calculated as 

« (X,) fdrj fdri'W{r,,7j,)gir^')(S{r^,eMv',0b)), (26) 



where qe{0b) = + T{lb ' '^)'^ ^\th ^f, « {Obxi^by, l)j and the functions g and W are given by equations (^) 
and (|2C|), respectively. In the second line on the r.h.s we have replaced the ionization fraction of X,., which 
generally depends on conformal time rj, with the average fraction, (Xg), between 77ion < rj < rjo assuming 
that Xa{vi) is a slowly varying function with respect to conformal time after a given reionized epoch Zjon- 
From equation (^, the ensemble average {5{6a,r])qe) can be rewritten in terms of third order moments of 
perturbations as 

{6{0a)qe{eb)) « 2(<5(l)(0„)[-^;f)^;f)+7^;(l)^;(^)](0b)) 

+ {5^^\ea)[-v'^^v'^^ +lv^^\^^\eb)) +0{{5^''>f). (27) 

Here we have used (^^^^ij^^^w'^^^) = and, from the fact that the propagation direction of CMB photons is 
almost parallel to z-direction in our setting, we have used 7^ • ~ , where Vz denotes the z component of 
bulk velocity field. Therefore, using equations (|l^) and (^, we can obtain 

X [4G'(-fc, k + k') + F{-k', k + k')] P{k)P{\k + fc'|)e'fc-L-('-^'»-'-'^^)-«fc^('?-'?')^ (28) 

where fc_L are (x, 2/)-components of k that are perpendicular to the line of sight, and r and r' denote the 
comoving angular distances from the present time to 77 and 77', respectively. One can readily see that the 
integration of the term e~'^^'^'^^'^ ^ over on the r.h.s of equation ( p8|) leads to the Dirac delta function 
5{'q—rj') as long as we are concerned with wavelength modes much smaller than the comoving Hubble distance 
characterized by 77. This leads to the consequence that dominant contribution to the cross-correlation function 
arises from structures of dark matter and ionized medium that are distributed at the same redshift and are 
separated by 6 on the observed sky. This is the so-called the Fourier-space analogue of Limber's equation 



(Kaiser 1992; Jaffe, & Kamionkowski 199? ). It is known that this fiat sky approximation causes errors no 



greater than 0(0.1%) at I > 1000 ( |Hu 2000b| ; [laffe, fc Kamionkowski 1998] ). Under this approximation, the 



cross-correlation function C"" {9) can be calculated as 



')0 



fk I dk I 

C^^\e)^{X,)l dr^g{7^)W{rj,f^.)D'D^ I^^S{kMk±r0), (29) 



»?ion 



with 



sik) . 2/ dk' rd^fV'^'^'-'^-''''''''^'^ 



JO 



A:2 + fc'2 + 2kk' cos 0a/1 - /x^ 

X [4G(fc, k', 0, /i) + F(fc, k', (j), fi)] P{k')P{\j F + fc'2 + 2kk' cos </> ^1 - A^^), (30) 

where we have assumed that k±^ is along the x-coordinate from the statistical isotropy and then fej^ • k' = 
k±k' cos(?!>cos0. 
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Hence, from equations (|2J) and (^9|) wc can finally derive the angular power spectrum of the cross- 
correlation between the convergence field and the secondary temperature fluctuation field by the quadratic 
Doppler effect of 0{v'^): 

Cf = (X,) iy^g{rj)W{rj, ^^)r-\^)b^ s{k = . (31) 

This is the first main result of this paper. It should be noted that the integrand quantity in C™'^ depends 
on the cosmological parameters and on the dark matter power spectrum. Namely, it is not affected by the 
unknown biasing relation between the dark matter and ionized medium distributions. In this sense, we can 
accurately predict the magnitude of the integrand quantity once we fix the dark matter power spectrum such 
as the CDM model and the cosmological parameters, which will be well constrained in the near future by 
other observations such as measurements of the CMB anisotropics, the weak Icnsing survey and so on. We 

2 

therefore expect that measurements of Ci " can set precise constraints on the average ionization fraction 
(Xe) between the reionization epoch and present. 



4.4. Density-Modulated Quadratic Doppler Effect of 0{v'^6f,) 

Next, we consider the secondary temperature fluctuation field induced by the density-modulated quadratic 
Doppler effect of A ^ 0{5eV^). From equation ([l7| ) we have 

h'^\e) - jdT^g(T^)X,5, [-v^ + 7(7 • vf] . (32) 

In the following, to relate 5e to the dark matter fluctuation field (5, we introduce the linear biasing parameter 
h expressed as 5e{x, rj) = b5{x, j]) for simplicity. Note that b generally depends on time and spatial position. 

It is relatively straightforward to calculate the cross-correlation function between the convergence field 

c 2 

and the temperature fluctuations field A : 



{X,b) jdv I dr,'W{v,r^.)g{rj'){SS){~v^ + 7{jf, ■ vf) 

4 



{X,b) Jdv jdr,'W{r^,r^Mmv,ea)5{r,',e,))^{v^), (33) 

where (X^b) is the average quantity of Xeb between the reionization epoch and present. After the similar 
procedure as used in the derivation of equation (^T[), we can obtain the angular power spectrum 

C^''^' « MX^b) jyvW{v,v*)9iv)r-'{v}D'p(^k = ^) v^^Jij), (34) 

where Vrms is the rms of the present-day bulk velocity field and can be expressed in terms of the linear density 
power spectrum as 

n , , D^(ri) fdk' , 
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Here the factor 1 /3 is due to the three-dimensional freedom of bulk velocity field such that Wrms represents 
the rms of one-dimensional motion. The most important result of equation is that the angular power 
spectrum C^' " is proportional to not quadratic but linear term of {Xpb) for all I and therefore could be 
both positive and negative, depending on the sign of {XJj) or more specifically on the sing of (6), since 
the integrand quantity in equation (jsj) is always positive. The quantity {XJ)) should be a crucial clue 
to understanding of how the ionized medium distribution is related to the dark matter distribution in the 
large-scale structure, which still remains large uncertainties both from theoretical and observational aspects. 
The measurements of Ci'' " can thus be a new direct probe of revealing the projected biasing relation. It 
is worth noting that the angular power spectrum of the OV effect depends on the quadratic term 0{X^b^) 
and is always positive irrespective of the sing of (Xeb). 

The simple assumption that the ionized medium traces the dark matter distribution allows us to place 
the upper limit on ^ using the nonlinear dark matter power spectrum, which has been well studied with 



N-body simulations for the full range of adiabatic CDM models. Equation (34) can be then rewritten as 



C^''^'' = iiXeb) rd7jWirj,v.)girj)r-'P^L (k=-,7j] <,(r?), (36) 



Vian 



r 



where P^h{k, rj) is the power spectrum of nonlinear density fluctuations for which we will employ the fitting 
formula developed by Peacock & Dodds (1996). Note that we here consider a possible case that electrons in 
reionized medium with non- linear density fluctuations are distributed along the line of sight moving according 
to a large-scale coherent bulk velocity field, as suggested by the CDM power spectrum. This is in fact an 
essential mechanism that the Ostriker-Vishniac effect can avoid cancellations between blue- and red shifts 
along the line of sight from the coupling of the density and velocity fields. Based on this consideration, in 
Figure ^ we use the non-linear power spectrum only for the calculation of C^'^^ and use the second-order 
perturbation theory for that of C^"" , because the latter arises from the ensemble average of the density 
fluctuation field and the bulk velocity filed on the mildly non-linear regime. 



RESULTS 



We present numerical results for the angular power spectra given by equations (|3|) and (|3j) for the 
ACDM model. We adopt a set of cosmological parameters of fimo = 0.3, flxo — 0.7, ilbo 0.05, h — 0.7, 



and (Tg = 1, which is consistent with the recent CMB experiments ( [Nctterfield, et al. 200l| ) and the cluster 
abundance (Eke, Cole, & Frenk 1996; Kitayama & Suto 1997). We simply consider the homogeneous 
hydrogen ionization model after a given reionization epoch in the range of 5 < Zjon < 30. In this case we 
will consider X,, = const, with respect to time during < z < Zion and equations (|^) and (|^) can be used 
for calculations of the optical depth and visibility function, respectively. Although more realistic model is 



2001), our model allows us to compute the 



a patchy or inhomogeneous reionization (e.g., Benson et al. 
magnitude of effects in the simplest way. 

We first estimate the order of magnitude of secondary temperature fluctuations by the quadratic Doppler 
effect. Using the ACDM model above, the present-day rms of bulk velocity field can be estimated from 
equation ( ^ ) as frms(^o) ~ 1-19 x 10^'^ in the unit of c = 1. Equation ([l7| ) tells us that the magnitude of 
temperature fiuctuations induced by the quadratic Doppler effect can be roughly estimated as A" ~ ~ 
10~^ — 10~^ for Zion = 10, since we have r ~ 10"'^ — 10~^ for the model considered in this paper. Similarly, the 
magnitude of temperature fluctuations induced by the OV effect can be estimated as A*-*^ ~ rVrmsS ~ 10~^, 
which is consistent with the numerical results done by Hu & White (1996) (see also Efstathiou 1988 and 
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Fig. 1. — Angular power spectra of the cross-correlation between the weak lensing field k and the secondary 
CMB fluctuations field A induced by the quadratic Doppler effects of ionized medium. These spectra are 
calculated under the simple assumptions that the gas traces the dark matter and the full hydrogen ionization 
after a given reionization epoch; X^, = 1 and h = 1 for the ionization fraction and the biasing parameter. 
The solid and dashed lines show contributions from the pure quadratic Doppler effect of O(u^) and the 

and (^J), respectively. We here considered five 
cases of the reionization epoch as Zion = 5, 7, 10, 15, 20. 



density- modulated effect of 0{v^5e) given by equations (31 



Jaffe, & Kamionkowski 1998). We thus conclude A" <C A^^ in the observed CMB map itself, and the 
magnitude of A^ is much below sensitivities of the satellite missions MAP and Planck Surveyor. 

Figure ^ shows the angular power spectra of the cross-correlations between the weak lensing field and 
the secondary temperature fluctuations induced by the quadratic Doppler effect of 0{v^) (dashed lines) 
and the density-modulated effect of 0{5ev'^) (solid), which are computed using equations ( ^T]) and (jsj), 
respectively. We here assumed h = \ for the unknown biasing parameter between the dark matter and 
electron distributions in the large-scale structure. Both the magnitudes of the Doppler effects of 0{v^) and 
0{5eV^) are comparable on / ^ 10^ in the case of no biasing 6=1. This is because leading contribution to 
the cross-correlation between the weak lensing field of 0{5) and the Doppler effect of 0{v^) comes from the 
fourth-order moments of linear perturbations like 0(5(^)^u*^^)^) resulting from the vanishing of the third-order 
moments = 0, while the cross-correlation from the density-modulated Doppler effect of 0{6eV^) 

also arises from the moments of (J*^^)?;^^)^). It is clear that the contribution from the density-modulated 
effect of Oiv^Se) dominates the effect of O(u^) on angular scales of I ^ 100. The total magnitudes of the 
angular power spectra peak Cf^^Z^/27r ~ lO^^*' around I w 1500 for the ACDM model and for the five 
cases of reionization epoch Zion = 5, 7, 10, 15, 20. The peak location in Z-space results from the shape of 
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CDM power spectrum. On the other hand, Cooray (2000) investigated the angular power spectrum of the 
cross-correlation between the weak lensing field and the thermal SZ temperature fluctuation field based 



on the halo approach of dark matter clustering (Komatsu, fc Kitayama 1999) and the simple model of 



intracluster gas distribution. The magnitude of the angular power spectrum from the SZ effect amounts to 
(jK,BZj2 ^ io~* at / « 2000 for the similar models as considered in this paper. The temperature fluctuation 
induced by the quadratic Doppler effect is over an order of magnitude smaller than that by the SZ effect 
because of A" ~ 10"^ - 10"^ and A^^ - 10"**. From this fact, it seems unlikely that the magnitude of 
the cross-correlation between k, and the quadratic Doppler effect is smaller only by ~ 10^^ than that from 
the SZ effect. This is because the cross-correlation with the Doppler effect comes from wide-range redshift 
structures such as < z < Zion, while the SZ effect has dominant contributions from low redshift massive 
dark halos such as clusters of galaxies with M ^ 10^'^ Mq at z ^ 0.5 ( Cooray 2000 ). These results therefore 



imply that, if the weak lensing fleld is a priori reconstructed by measurements of lensing effects on the 
CMB, the cross-correlation could improve the S/N of the quadratic Doppler effect, even though the angular 
power spectrum itself of the Doppler effect has a small power of CP^P ^ 10~^^ only. This is an analogous 
statement to the cross-correlation between the primary CMB temperature fluctuation and polarization fields 



(e.g., |Hu, fc White 199?! ) 



To explicitly illustrate how the cross-correlation considered could be useful to reveal the biasing relation 
between the dark matter and electron distributions, we here employ a toy model of the antibiasing relation, 
motivated by the reionization scenario of low-density region prior to the high-density region as suggested by 
MHR. The homogeneous full hydrogen ionization with X^, — 1 seems to lead to 6^ — 5b from Uf, — nn , which 
seems inconsistent with the antibiasing relation, and thus we here assume the partly homogeneous ionization 
fraction of X,, = 0.5 after Zjon- Although the antibiasing process depends on how reionization proceeds in 
the large-scale structure in detail, we consider the simplest model of linear and constant antibiasing relation 
given by 6e — bS with b = —0.8. It should be then noted that the condition of 6e > —1 leads to S < l/\b\. 
Accordingly, we have taken into account the range of 6 for the calculation of the ensemble average of {6eS) 
used in equation (^4|). Figure ^ shows the result of this toy model, and one can readily see that we should 
observe the negative total angular power spectra (solid lines) of the cross-correlation between the weak lensing 
field and the secondary temperature fiuctuation field induced by the Doppler effect on I ^ 2000. Note that 
this critical scale of I depends on the value of the antibiasing parameter 6, because the cancellation between 
the spectra from the quadratic Doppler effect of 0{v'^) and the density-modulated effect of 0{v'^Se) also 
depends on b. As expected, the cross-correlation function could thus be a direct probe of the antibiasing 
relation, which cannot be attained by measurements of the angular power spectrum of the 0{V) Doppler 
effect such as the OV effect because of the quadratic dependence on the biasing relation. 

Figure ^ shows the maximal estimate of the angular power spectrum of the cross-correlation using 
the nonlinear dark matter power spectrum under the simplest assumption that the gas traces dark matter 
distribution in the large-scale structure. We here again considered the simplest model of X^ = 1 and 5 = 1 as 
in Figure |l|. The nonlinear enhancement could be important on small angular scales of I }t 10'^. However, it 
is not still clear if the electron distribution traces dark matter distribution in the nonlinear regime, because 
physical processes of the gas pressure and radiative dynamics should also play an important role in such 
nonlinear regions (e.g., Gnedin 2000). For this reason, the curves in Figure |^ provide us with upper limits 



on a realistic angular power spectrum of the cross-correlation. 
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Fig. 2. — A toy model of cross-correlation angular power spectrum for the simplest model of antibiasing 
relation between the dark matter and electron distributions given by 6 = —0.8. We here considered the 
reionization epochs of Zion = 5, 7, 10, 15, 20 and assumed X,, — 0.5 for the ionization fraction after the 
respective reionization epoch. The dashed and dotted-dashed lines show the angular power spectra of cross- 
correlation between the weak lensing field and the temperature fluctuation field induced by the Doppler 
effects of O(u^) and 0{v'^5e), respectively, as shown in Figure [I[ Solid lines show the total angular power 
spectra that we would actually observe, and the spectra become negative on I ^ 2000 for this model. 



6. DISCUSSION AND CONCLUSIONS 



In this paper, we have investigated the cross-correlation between the weak lensing field and the secondary 
CMB temperature fluctuations induced by the Doppler effect of the ionized medium in the reionized epoch 
of the universe. The main result of this paper is to develop a formalism for calculating the cross-correlation 
function and its angular power spectrum based on the small angle approximation. The leading order to the 
cross-correlation comes from the secondary CMB anisotropics due to the quadratic Doppler effects of the bulk 
motions, since the cross-correlation with the 0{v) Doppler effect such as the OV effect depends linearly on 
the bulk velocity field and thus suffers from cancellations between positive and negative contributions from 
the bulk velocity field for the ensemble average. There are two main contributions to the cross-correlation 
function: one is contribution from the pure quadratic Doppler effect of 0{v^) and the other is that from the 
density- modulated effect of 0(u^(5e), where 5e is the electron density fluctuation field. For the CDM models, 
the latter contribution dominates the former on angular scales of ^ 100 if the ionized medium traces the 
dark matter distribution. We found that the cross-correlation of 0{XeV^5) from the Doppler effect of 0{v^) 
is sensitive to the average ionization fraction {X^) between the reionization epoch and present in the sense 
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2x10 




Fig. 3. — Maximal nonlinear enhancement of the angular power spectra of the cross-correlations for the 
same models {X^ = 6 = 1) as shown in Figure |l]. Solid and dashed lines demonstrate the total angular 
power spectra taking into using the linear and nonlinear dark matter power spectrum, respectively. These 
power spectra includes both contributions from the cross-correlations of the weak lensing field with the pure 
quadratic Doppler effect of 0{v'^) and the density-modulated effect of 0{v'^5e) as shown in Figure |l|. We 
used equation ( p6| ) for the calculation of the nonlinear density-modulated spectrum. Nonlinear effect could 
be important at I ^ 1000. 



that the bulk velocity and the weak lensing field depends only on the dark matter fluctuations, through the 
gravitational instability scenario, and on the cosmological parameters, which will be well constrained by the 
future precise data such as the primary CMB anisotropics and the weak lensing survey. On the other hand, 
we showed that the cross-correlation of 0{Xev'^6eS) between the weak lensing field and the density-modulated 
quadratic Doppler effect of 0{v'^Se) depends linearly on the biasing relation between the dark matter and 
free electron distributions, and therefore its angular power spectrum could be both positive and negative, 
depending on the biasing and antibiasing, respectively. If the low density region of IGM is ionized prior to 
the high density region in the reionization history as suggested by MHR (see also Gncdin 200C| ), we would 
measure the negative angular spectrum of the cross-correlation. Our results thus offer a new opportunity 
to understand the unknown problem how reionization proceeds in the inhomogeneous intergalactic medium 
connected to the dark matter clustering in principle, although the magnitude of the signal is very small. 

Let us briefly comment on other sources of cross-correlation between the weak lensing field and the 
observed CMB map. One significant source is the thermal SZ effect, and it has been indeed shown that the 
cross-correlation could have the magnitude of (A^^k) ^ 10^^ for the similar CDM model ( Cooray 2GGC ) as 
investigated in this paper. The SZ effect thus produces much larger contribution than the cross-correlation 



with the Doppler effect estimated as (A 



DP. 



10 



-10 



However, since the thermal SZ signal arises mainly 
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from hot plasma in clusters of galaxies at low redshift (z ^ 0.5), it is expected that we can remove this signal 
from the cross-correlation maps, combined with the X-ray data or the SZ maps. The other possible sources 
are the CMB foregrounds such as synchrotron radiation from extragalactic dust that may correlate with the 
dark matter distribution in the large-scale structure. Isolating the signal shown in this paper observationally 



from the foregrounds will be a great challenge by using the different frequency-dependence (Tegmark, et al 
2000] ), though this problem is beyond the scope of this paper. 
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